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FLUTTER CALCULATIONS IN THREE DEGREES OF FREEDOM 

By Tiikodore The odors en* and I. E. Garrick 



SUMMARY 

The present paper is a continuation of the general study 
of flutter published in AVICI Reports Xos. 496 and 6S5. 
The paper is mainly de retted to flutter in three degrees of 
freedom (bending, torsion, and aileron), for which a 
number of selected cases have been calculated and pre- 
sented in graphical form. The results are analyzed and 
discussed with regard to the effects of structural damping, 
of fractional-span ailerons, and of mass-balancing. Ute 
analysis shows that more emphasis should be put on the 
effect of structural damping and less on mass-balancing. 
The conclusion is drawn that a definite minimum ambunt 
of structural damping, which is usually found to be 
present, is essential in the calculations for an adequate 
description of the flutter case. TJieoretieal flutter pre- 
dictions are thus brought into closer agreement with the 
facts of experience. 

A brief discussion is included of a particular biplane 
that had experience d flutter at about 200 miles per hour. 
Some simplifications hare been achieved in the method of 
calculation. 

INTRODUCTION 

Since the publication of the previous flutter papers, 
the necessity of considering complete cases of three 
degrees of freedom including the effect of structural 
damping has become evident. The purpose of the pres- 
ent paper is therefore to present such extensions of gen- 
eral applicability. The calculations herein reported are 
directly based on methods already given in references 1 
and 2. The earlier papers deal, to some extent, with 
cases of three degrees of freedom and also indicate that 
the internal structural damping in some cases has a 
great effect on the flutter velocity; a small value of the 
internal damping may suffice to bring the flutter veloc- 
ity from nearly zero to a normal value. Thus, in order 
to obtain better agreement with practice, the existence 
of a certain amount of internal damping must be 
recognized. 

A separate investigation on the subject of hysteresis 
in airplane structures, which has been conducted in the 
meantime and will be report-oil in detail elsewhere, 
shows that a significant amount of internal damping 
(g a > 0.01) is present, usually with considerable mar- 
gin. This low value of g a ~ 0.01 is found to be effec- 



tive in smoothing out the low-velocity flutter values 
appearing in flutter curves calculated for the case of 
zero internal damping. A similar effect of different 
origin is the so-called fractional aileron-span effect. 
This effect was noted in reference 1 for binary cases and 
is here also treated for ternary cases. Strangely 
enough, a reduction in the length of the aileron from 
that of the full span to a shorter length has a dispropor- 
tionally large effect on the flutter velocity. Thus, the 
calculated flutter speed for a full-span aileron may be of 
a low value; whereas, for a half-span or even a three- 
quarter-span aileron, it may be nearly normal. 

It is of interest to note in connection with the study 
of three degrees of freedom that the addition of the 
third degree is the cause of a reduction in the flutter 
speed based on only two degrees. If a control surface is 
mass-balance* I, is reasonably stiff, and a certain mini- 
mum amount of torsional damping is present, the 
bonding-torsion value of the flutter speed will be closely 
approached. 

The following ^tudy originated in an investigation of 
a certain biplane in which flutter had been experienced 
on a number of occasions. Two of these biplanes were 
made available at Langley Field for the purpose of the 
investigation. These biplanes were subjected to the 
conventional vibration tests in order to obtain the 
flutter parameters, and the flutter speed was calculated. 
These calculations were used as the nucleus in the fol- 
lowing study of flutter in three degrees of freedom. For 
readers particularly interested in the biplane mentioned, 
an appendix (appendix C) has been prepared. 

It should further be mentioned that some simplifica- 
tion has been achieved in the method of calculation. 
This simplification is based on an analogy with Sylves- 
ter's method of elimination and reduces quite noticeably 
the labor of calculating the flutter speed for three 
degrees of freedom. Appendix B presents a summary of 
this method. 

RESULTS 

The results of the flutter calculations are presented 
in figures 1 to 40. In tables I to IX the constant 
parameters and the variable parameters arc arranged 
to serve as a key to the figures. In order further to 
assist the reader in the stutly of the curves, a brief 
description of the figures will be given. 
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It will be noticed that the ordinate for all the curves 
is the flutter speed in the coefficient form v/b<a a . The 
product J>u> a is thus used as a reference velocity through- 




i.6 /.a ao 



Figure 1.— Flutter coefficient r/&«« against frequency ratio wj/w. for several values 
of the aileron unbalance, z fi . j.,0.2; r«», 1; no damping. 

out. The symbols used in this paper are defined in 
appendix A. 

The figures are arranged according to the values of 
r a 2 : figures 1 to 11, r a 2 ~ 1 (biplane case); figures 12 to 
28, rv = 0.5; figures 29 to 36, r o 2 =0.25 (monoplane 
case). Within each group a further arrangement is 



which corresponds to the bending- torsion binary 
flutter value. 

Figure 2 shows the effect of the torsional structural 
damping coefficient 0 O =O.O1 on some of the curves of 
figure 1. Note that the dip in the flutter curves is now 
eliminated and that the flutter coefficient does not 
differ by much from its bending- torsion value. 

Figure 3 shows the individual effects of the structural 
damping coefficients g a , g fi , and g h on the flutter co- 
efficient for the constant parameters 0.002 and 
Wwa = 0.833. Note that g a has the greatest effect in 
increasing the flutter speed. 

The parameters for the next set of curves (fig. 4) 
differ from thoso of figure 1 only in the value of x«, 



.4 




Asymptotes 



0 

(a) -0.004. 

(b) -0.002. (c) ij- -0.002. 

Figure 2.— Flutter coefficient r/6oj a against frequency ratio vf!<*» a with and without structural torsional damping. z -f 0.2; r«* t 1. 



i.Q 2.0 



made according to the value of k, the wing-density 
parameter. 

Figure 1 shows a number of curves plotted against 
the aileron frequency ratio o)?/o) a) with « a thus used as a 
reference frequency. The wing bending-frequency 
ratio o? A /^a is kept constant. The curves differ only 
in the value of j% which determines the degree of aileron 
mass balance. Note the low dips present nearw^/w a = 1.0 
and the shifting of these low spots with the value of 
x$. .All the curves approach an asymptote for w^/co a -> oo t 




Figure 3.— Flutter coefficient apainst structural damping coefficients Sfi, 
and ga. «0/u«, 0.833; if, 0.002; r«>, 1. 



which is now 0; that is, the center of gravity of the main 
wing coincides with the elastic axis at the 40-percent - 
chord position. Again, for values of x& of 0 and 0.002, 
low dips exist near w^/w a = 1.0. For 0.002, the 

low dip does not exist. The bending-torsion flutter 
value at w$jw a —v> is considerably increased over that 
for x a =0.2 in figure 1. 
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Asymptotes 
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0 .3 .4 .6 .8 LQ !.? 1.4 i.6 i.8 e.Q 

Figure 4.— Flutter coefficient p/tai a against frequency ratio w*/a> a for several values 
of the aileron unbalance ij. 0; r«», 1; no damping. 
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8 1.0 1.2 1.4 1.6 l.t 2.0 



(a) it-o, ** m -0.002. 

Figure 5. -Flutter coefficient r/6« a against frequency ratio with and without structural torsional damping. j«, 0; r„», 1. 

L6 




0 .2 .4 .6 .8 I.Q 1.2 1.4 1.6 1.8 2.0 

Figure 6.— Flutter ccK^fficitnt r;bu a against frequency ratio wafv a for several values 
of tbe wing unbalance x*. with and without structural torsional damping, it, 
0.002; r a », 1. 
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0 .2 .4 .6 .8 I.Q 1.2 1.4 1.6 I.Q 2.0 

Figure 7.— Flutter coefficient r/6a> a against frequency ratio with and without 
structural torsional damping. r«, —0-1; is, -0.005; r«», 1. 
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Figure 8.— Flutter coefficient r/frw« against frequency ratio u}/w a showing the effect 
of partial-span aileron coefficient {, is, 0.002; r«', 1. 
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Figure 9.— Flutter coefficient vjbui* against partial-span aileron coefficient u>a/u>«, 
0.833; r«», 1. 




Figure 10.— Flutter coefficient r/6w„ against frequency ratio uj/u,, showing the com- 
bined effect of structural damping coefficient g m and partial-span aileron coeffi- 
cient I. Xfi, 0.002; r«», L 
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Figure 1 1.— Flutter coefficient r/bu> a against frequency ratio u«/« a showing the com- 
bined effect of the structural damping coefficient g a aDd partial-span aileron coeffi- 
cient i. is, 0.002; r a *. 1. 
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Figure 5 shows the effort of the torsional structural 
damping coefficient <7 a = 0.01 in increasing the value of 
the flutter speed. Figure G gives several curves for a 




Fig i re 12.— Flutter coefficient r/ta« against frequency ratio ug/u, for several values 
of the aileron unbalance ifi. i« ( 0.2; r„i, 0.5; no damping. 



curves presented in figure 4 (x« = 0, x$ — 0.002). The 
effect of <7 tt = 0.01 is shown for comparison. It is inter- 
esting to observe that in the range wj/<u a <Cl.O the effect 
of £ is significant. In the comparison of this case with 
figure 10 (x a = 0.2), it appears that £ = 0.8 is of more 
influence on the case ar a =0 while # a =0.01 is more 
effective on the case z a =0.2. 

The next set of figures (figs. 12 to 28) has been cal- 
cidatcd with r a 2 =0.5. Figure 12 is similar to figure 1 
and shows the flutter-speed coefficient plotted against 
aileron frequency ratio for several values of xp. The 
effect of structural damping is included in figure 13. 
Figure 14 is a cross plot (similar to fig. 3) against the 
structural damping coefficients g aj g$, and g h . Figure 
15 extends the cases given in figures 13 (a) and 13 (c) to 




:.0 Id i.4 id t.S so 



(a) j tf = 0.002. 
(C) = -0.002. 



(b) r*-0. 

(d) x*- -0.005. 



Figure 13.— Flutter coefficient r.'b^ against frequency ratio u^/w. with and without structural torsional damping, r., 0.2; r a \ 0.5; u 0.607. 



constant value of x$ of 0.002 and for different values of 
x a (0.2, 0, and —0.1), with and without structural 
damping. 

Figure 7 represents a case for which x a =— 0.1 and 
xp= — 0.005. Case 1 (bending-torsion) is completely 
stable- 
Figure S shows the effect of £, the partial-span aileron 
coefficient. The curve $=1.0 is taken from figure 1 
(i0 = 0.002) and is the case of the full-span aileron. 
Note that even a small reduction to J = 0.8 has a marked 
favorable effect, especially in the range of frequencies 
1.0. As £-»0 (no aileron), the curves approach 
the bending-torsion flutter value. 

Figure 9 represents a plot against £ for a constant 
value of Wjs/wa of 0.833. Figure 10 is intended to show 
a combined effect of £ = 0.S and <7«=0.01. For com- 
parison the separate combinations £ = 1.0, g a — 0; 
£=1.0, g a — 0.01; and £ = 0.8 ? g a ~0 are also shown. 
Figure 11 shows the effect of £=0.8 on one of the 



include other values of the frequency ratio w A /u> a . Fig- 
ure 16 represents a case of a lighter wing for which < is 
0.25 instead of 0.2. The value of x & is 0.002; curves 
with and without structural damping are given. Figure 
17 has the same conditions presented in figure 16 except 
that x a is equal to 0 instead of 0.2. 




(a) ufilu at 0.6O7. 



(b) «,/« fl , 0.S33. 



Figure 14.— Flutter coefficient r/6a> a against structural damping coefficients g a , g$ t 
and Qk. ifl, 0.002; r*\ 0.5. 
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KO L2 1.4 1.6 1.8 2.0 



(a) xs, 0.002; 0.316. 
(c) I,, -0.002; wi/w, 0.316. 

Figure 15.— Flutter coefficient vfbw. against frequency ratio with and without structural torsional damping. t«, 0.2; r«», 0.5; (cf. figs. 13 (a) and 13 (c)). 



(b) Xf, 0.002; on/—, 1.0. 
(d) x fit -0.002; 1.0. 
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Figure 16 —Flutter coefficient tfbw a against frequency ratio with and without 
structural torsional damping. x«, 0.2; Xfi, 0.002; *, 0.25; r.», 0.5. 
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0 .<? .6" .8 t.O 1,2 1.4 1.6 1.8 2.0 

Figure 18.— Flutter coefficient v(bu> a against wing bending-frequency ratio 
u fi fu a , 0.5; if, 0.002; «, 0-25; r«>, 0.5; no damping. 




.8 t.O 1.2 14 1.6 1.8 2.0 

Figcbe 20.— Flutter coefficient vfbua against wing bending-frequency ratio for 
several values of the aileron unbalance 19. upfe/a, 0; «, 0.2; r«>, 0.5; no damping. 
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0 .2 .4 .6 .8 (.0 1.2 1.4 1.6 1.8 2.0 

Figure 17.— Flutter coefficient p/6«« against frequency ratio with and without 
structural torsional damping. x m , 0; if, 0.002; 0.25; fa 1 , 0.5. 




Figurk 19.— Flutter coefficient p/6w« against wing bending-frequency ratio 
wj/wa, 1.0; xf, 0.002; *, 0.25; r«>, 0.5; no damping. 
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0 .2 .4 .6 .8 1,0 1.2 1.4 1.6 l.d 2.0 j 0 .2 A .6 .8 1.0 1.2 IA f.6 1.8 2.0 

Figure 21.— Flutter coefficient against'frequency ratio <##/«• with and without Figure 22.— Flutter coefficient vfb» m against frequency ratio with and without 
structural torsional damping, zp, 0.002; *, 0.125; r«>, 0.5; 0.807. structural torsional damping, x$, -0.002; », 0.125; r«», 0.5; 0.507. 





0 005 . 01 .0/5 . 02 .025 . 03 .035 .04 



Figure 24.— Flutter coefficient vjbu a against structural damping coefficient g». 
0.316; x», 0.002; «, 0.125; r«», 0.5. 
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Figure 25.— Flutter coefficient v(bw 9 against frequency ratio u^/u. showing the com- 
bined effect of the structural damping coefficient g, and partial-span aileron coeffi- 
cient <. ip, 0.002; r, 0.083; r J, 0.5. 



FLUTTER CALCULATIONS IN THREE DEGREES OF FREEDOM 




(a) «*/««,, 0.316. ( b ) «4/„. t 0.607. 

(c) 1.0. 

Figure 26. -Flutter coefficient v/b» m against frequency ratio o*/«« with and without structural torsional damping. x ft 0; r«, 0.2; *, 0.083; r«» ( 0.5. 
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.8 io /.<? i4 t.6 /.e 2.0 

Fiou*x 28.— Flutter coefficient tfbw. against frequency ratio 1.0; r. f 0; 

r, 0.26; r m \ 0.5; no damping. 
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.8 t.O 1.2 IA 1.6 f.8 2.0 

Figcbi[29.— Flutter coefficient r/6«. against frequency ratio for two values of 
aileron unbalance ij. r„», 0.25; no damping. 




(a) Xfi, 0.002. (b) z$, -0.002. 

FiGXTEE 30.— Flutter coefficient r/6w« against frequency ratio with and without 
structural torsional damping. u\ 0.25; 0.007. 
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Figure 31.— Flutter coefficient t?/6^. against frequency ratio for two values of 
the partial-span aileron coefficient f. Xfi, 0.002; r«', 0.25. 




(9 .i> .6 .8 1.0 L2 1.4 /.€ 1.8 O .2 A .6 .8 t.O L2 L4 1.6 1.8 2.0 

v 0/ 'u> m 

(a) Xfi, 0.002; wa/««, 0.318. (b) X$, 0.002; 1.0, 
(c) xfi, -0.002; 0.316. (d) Xfi, -0.002; wW« a . 1.0. 



Figuee 32.-Flutter coefficient p/6o»« against frequency ratio with and without structural torsional damping. (See also fig. 30.) «, 0.2; r.i, 
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Figure 18 is a plot of the flutter coefficient against the 
wing bending-frcqueney ratio, for a constant value of 
cififwa = 0.5. The case w A /w fl = oo corresponds now to 
the binary case, torsion-aileron. The branch repre- 
senting essentially this case is easily evident. Figure 
19 differs from figure 18 only in the value of «j/w a , 
which is now 1.0. The branch representing torsion- 
aileron is now gone. (The small singular branch on the 
axis near u A /w„ = 1.1 can be shown to disappear com- 
pletely with a very small amount of friction.) 

Figure 20 differs from figure 18 in the value of 
which is now 0.2, and also in the value of which 
is now 0. In addition, several values of x$ have been 
employed. Note that the aileron-torsion branch 
beyond 1.0 exists only for the largest unbal- 

ance, T^O.002. 



Figure 21 differs from the parallel cases shown by 
curves 7$= 0.002 in figures 12 and 10 only in the value 
of which is now 0,125; that is, it represents a heavier 
wing or a higher altitude. Note that xj=0.002 does 
not eliminate the torsion-aileron branch. The effect 
of 0 a = O.Ol produces a flutter curve, the ordinate of 
which is remarkably near the bending- torsion value. 

Figure 22 differs from figure 21 in the value of z$, 
which is now —0.002. The low dip near ^/a? a = 1.0 
is eliminated for a value of g a — 0$\. Figure 23 
extends the cases of figures 21 and 22 to two other 
values of the frequency ratio w A /w a . 

Figure 24 is a plot of the flutter coefficient against 
g a for the constant value of j^ = 0.002 and w^/w a = 0.316. 
(See fig. 21.) Note that the torsion-aileron branch 
is gradually eliminated and vanishes for g a ^ 0.006. 




(a) ut/ua, 0.318. 0>) 0.007. 

(C) <•>»/<•»«, 1.0. 

Fiqube 33.— Flutter coefficient tjbu a against frequency ratio wpjw* with and without structural torsional damping. «, 0.125; r«», 0.25J ifl, 0.002. 




.2 A .6 .8 1.0 1.2 i.4 L6 i.8 20 

(a) 0.315. (b) 0.607. 

(C) OlkJUa, 1.0. 

Figure 34.— Flutter coefficient f/&«« against frequency ratio with and without structural torsional damping. «, 0.125; r«», 0.25; if, -O.0O2. 
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(a) «*/»., 0.316. 
(c) 0.707. 



(b) 0.5. 
(d) 1.0. 



Fiours 35 —Flutter coefficient p/6«« against frequency ratio with and without structural torsional damping. *, 0.0752; r«*, 0.25; xj,Q» 



Figure 25 represents a still heavier wing (k = 0.083). 
This curve shows that j^ = 0.002 does not eliminate 
either the torsion-aileron branch or the bending- 
aileron branch for low values of w^/to a . The value 




(a) K t 0.25. 

(b) *, 0.1. 

Figfei 36— Flutter coefficient r/6«« against frequency ratio a>j/W r«, 0.2; r«», 0.25. 



{—0.7 as shown eliminates the low branches. The 
value <j a =0.02 eliminates the torsion-aileron branch 
but has little influence on the bending-aileron branch. 
Figures 26 and 27 represent similar cases with x& = 0 
and with several values of the frequency ratio o> h /a> a . 
In the cases represented by figure 26 (r a =0.2) the 
center of gravity of the wing is at 50-percent chord 
and for those of figure 27 (x a =0) the center of gravity 
is at 40-percent chord. 

Figure 28 represents a case in which k = 0.25, x„=0, 
and o> /l /a> a =1.0. The figure shows that the bending- 
torsion flutter branch is eliminated and only the 
torsion-aileron branch exists. This branch can also 
be eliminated by increasing the value of g a . 




Figitke 37.— Flutter coefficient p/6« a against the structural damping coefficient g m In 
an tisym metrical cases for several values of the aileron unbalance ifi. 0; 

uk/u«i 0. t 
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Figure 38- Flutter coefficient r/^ a against partial-span aileron coefficient £ in the antisymmetrical cases for several values of the aileron unbalance if. 0; «»/««, 0. 



The next set of figures (figs. 29 to 36) have been 
calculated with r a 2 = 0.25 (monoplane case). Figure 
29 shows the flutter coefficient plotted against ^/co a 
for two values of 0.002 and -0.002. The effect 
of structural damping, ^ a = 0.01, is shown in figure 30 
and the effect of the partial-span aileron coefficient £ 
is shown in figure 31. Figure 32 extends the cases 
of figure 30 to other values of the bending-frequency 
ratio u h :cc a ; figures 33 and 34 represent parallel cases 
for a heavier wing, * = 0.125. 

Figure 35 represents a monoplane case with 
parameters based on a modern heavy pursuit airplane. 
For completeness, several curves are shown with dif- 




.004 



Figure 39.— Flutter coefficient r/bw* against aileron unbalance z$ in the antisymmet- 
rical cases for several values of the partial-span aileron coefficient {._u>0;W, 0; 



ferent values of the bending-frequency ratio 
Figure 36 is based on the parameters for a modern large 
airplane. Two values of * are presented: 0.25 and 0.1. 

The rest of the figures were • calculated for two 
constant values: w^/w a ^0 and oj ft /oj a =0 (antisym- 
metrical flutter cases). Figure 37 shows the flutter 
coefficient plotted against g a for four values of x & 
(0.004, 0.002, —0.002, and —0.006). It is observed 
that the effect of g a is quite significant. Figure 38 
shows the flutter coefficient plotted against £ for the 
same values of ^ that were used in figure 37. The 
effect of £ in figure 38 is rather large. Figure 39 is a 
cross plot of figure 38. with x$ as the abscissa. Figure 
40 is a plot of the flutter coefficient against g a for three 
values of rj (1, 0.5, and 0.25) and for two values of 
(0.002 and -0.002). 
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Figure 40.— Flutter coefficient r/6w a against the structural damping coefficient g a 
in the antisymmetrical cases for three values of r« J and two values of aileron un- 
balance ifi. 0; wk/wa, 0. 
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DISCUSSION 

The first noteworthy observation in the case of three 
degrees of freedom is the distinct dip in the flutter 
curve at values of u$/w« somewhat less than unity 
when structural damping is neglected. Apparently the 
aileron under these circumstances is very nearly in 
mechanical resonance with the wing in torsion. It is 
further observed that the flutter velocity remains 
rather low in this range of values of the aileron fre- 
quency. Since the aileron frequency in most practical 
cases is definitely less than that of the wing torsion, the 
region below r unity is of the most significance. 

There are two types of aileron response: One type 
corresponds to symmetrical wing motion and the other 
type corresponds to an an ti symmetrical motion. The 
frequency of the first type is of the order of one-half to 
three-fourths of the torsion frequency and the fre- 
quency of the second type is zero. It is noted that the 
elimination by mass-balancing of flutter resulting from 
the symmetrical type of response may be difficult, 
particularly if the aileron frequency is close to the wing- 
torsion frequency; whereas, the antisymmetrical type 
is more favorably affected by normal mass-balancing 
of the aileron. It is also to be noted that the wing 
damping is unusually effective in removing the dip in 
the flutter curve. Indeed, for comparatively light 
structures a value of the torsional damping coefficient 
g a of 0.01 brings the flutter velocity almost back to its 
full bending-torsion value. Significantly, the torsional 
damping seems to be the most effective. Heavier 
structures appear to be less susceptible to the effect of 
damping. In fact, a larger value of g a is needed and 
apparently it may be necessary also to provide damping 
in one or both of the other degrees of freedom (fig. 25). 

A partial-span aileron has a rather profound effect 
on the dip in the flutter curve, which is similar to the 
effect of the damping. A reduction of the effective 
aileron length £ from 1.0 to 0.8 practically restores the 
normal value of the flutter speed. 

It is rather evident from the present study that the 
effect of mass-balancing has been overemphasized in 
the earlier literature. Of significance is the fact that 
a pronounced dip exists in the flutter curve even for 
an overbalanced aileron (fig. 1). The aileron balancing 
seems to become most effective for the case in which 
the wing itself is overbalanced (fig. 7). This case is 



only of academic interest. Overbalancing alone docs 
not present a solution of the general case of three 
degrees of freedom; the appropriate value of the flutter 
speed cannot be obtained solely by any practicable 
method of balancing. 

On the other hand, the greatest beneficial effect of 
damping is obtained for the unbalanced, that is, the 
normal wing (fig. 6). Only in this case is the full 
bending-torsion value nearly reached. In the range of 
frequencies <o^/w a <C 1 the flutter speed of the overbal- 
anced wing remains much lower than that of the nor- 
mal wing. It is further noted that the beneficial effect 
of aileron balance is small when a small amount of 
damping is present (fig. 2). 

For the antisymmetrical case with no damping 
present, w A — 0, it is observed that the balancing of the 
aileron is more effective. For a given value of the 
torsional damping coefficient (^«=0.01) the gain from 
balancing is not large. The effect of the fractional 
aileron is very marked. At {=0.8 the flutter velocity 
equals the torsion-bending value independently of the 
balance coefficient. 

CONCLUSION 

It has been shown that mass-balancing is of less 
significance than has heretofore been attributed to it. 
The profound effect of internal structural damping has 
been shown. For the normal, unbalanced wing a small 
amount of damping removes the dip in the flutter curve 
and substantially yields the torsion-bending value of 
the flutter velocity. The large beneficial effect of the 
fractional-span aileron has been indicated. These state- 
ments apply to light, low-density structures and apply 
to a lesser degree as the wing density is increased. 
Because of the complexity of the problem, too general 
conclusions cannot be safely made and detailed calcu- 
lations of individual cases are still needed. The in- 
cluded graphs, which cover a fairly representative field, 
should be of value for specific studies and should furnish 
numerical solutions in a number of cases. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 7, 194L 



APPENDIX A 



LIST OF NOTATION 



a 

& 
h 
b 
a 



angle of attack (fig. 41) 

aileron angle (fig. 41) 

vertical distance (fig. 41) 

half chord, used as reference unit length 

coordinate of elastic axis (also called axis 
of rotation or torsional axis) (fig. 41). 
Location of elastic axis in percentage 
total chord measured from leading edge is 



100 



1+a 



- or a= 



2 (elastic axis) 

Too 



-1 




of en ti re wing ^ 



e.g. of aileron- 

Figure 41.— Half chord b is used as the unit length. The positive directions of 
a , fi, and h are indicated by arrows. Note that a is measured from midchord and 
z a is measured frum the elastic axis positive to the right. Also note that Zf is a 
♦'reduced' ' parameter and not the actual distance from the hinge to the center of 
gravity of the aileron. 

c coordinate of aileron hinge axis (fig. 41). 

Location of aileron hinge axis in percent- 
age total chord measured from leading 

edge is 



100 



1-fc 



or c- 



2 (aileron hinge) 
100 " 



-1 



P 

M 



irpb- 

= M 



x "~Mb 



mass of air per unit volume 
mass of wing per unit span length 
ratio of mass of cylinder of air of diameter 
equal to chord of wing to mass of wing, 
both taken for equal length along the 
span; this ratio may be expressed as k= 
0.24 (b 2 /W) (p/po) where W is weight in 
pounds per foot span, b is in feet, and 
p 'po is ratio of air density to standard air 
location of center of gravity of wing-aileron 
system measured jrom a (fig. 41); S a , 
static moment of wing-aileron per unit 
span length referred to a. Location of 
center of gravity in percentage total 
chord measured from the leading edge is 
l+a+r. 



100 



a+x a 



or 



2 (center of gravity) 
100 



-1 



x >~Mb 



C a 

c, 



- /S 



0)0 

'la 

w * \ A? 
t 

v 

. dib 



i/k 

vjbo) a 

9* 9? 9h 



reduced location of center of gravity of 
aileron referred to c (fig. 41); S$, static 
moment of aileron per unit span length 
referred to c. M refers to total wing 
mass and not to mass of aileron alone 

radius of gyration of wing aileron referred 
to a (fig. 41); Ia f moment of inertia of 
wing aileron about elastic axis per unit 
span length 

reduced radius of gyration of aileron re- 
ferred to c (fig. 41); moment of 
inertia of aileron about c per unit span 
length 

torsional stiffness of wing around a per 

unit span length 
torsional stiffness of aileron around c per 

unit span length 
stiffness of wing in bending per unit span 

length 

natural angular frequency of torsional 
vibrations around a in vacuum (co a = 
2t/«, where / a is in cycles per sec) 

natural angular frequency of torsional 
vibrations of aileron around c 

natural angular frequency of wing in 
bending 

time 

speed of forward motion 

angular frequency of wing vibrations 

reduced frequency = number of waves in 

wake in a distance equal to scmichord X 

2t 

reduced wave length —length of one wave 
of wake in terms of a distance equal to 
semichord X 2ir 

flutter-speed coefficient 

structural damping coefficients; rg cor- 
responds approximately to the usual 
logarithmic decrement 

partial-span aileron coefficient. Note that 
this coefficient is not the geometric 
ratio but an "effective" value of the 
order ot[ff(a)dz] 2 /ff 2 (a)dz, where the 
integral in the numerator is taken over 
the aileron span and that in the de- 
nominator is taken over the full span; 
f(a) represents the spanwise amplitude 
of (flutter) torsion mode 
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APPENDIX B 



METHOD OF ELIMINATION AS APPLIED TO FLUTTER CALCULATIONS 



The treatment of the flutter problem (references 1 
and 2) leads to the simultaneous solution of two equa- 
tions. The degree of each of these equations in the 
general case of three degrees of freedom (flexure, tor- 
sion, and aileron) is three. If, in addition, the effect of 
a tab motion or a float is desired, the degree of the 
equations may be more than three. The numerical 
calculations involving the plotting of roots becomes 
laborious and time-consuming. A method of elimina- 
tion for obtaining common roots of two simultaneous 
equations may be used, which does away with the 
necessity for any root extractions. (See, for example, 
reference 3.) The procedure results in the saving of 
considerable effort, particularly when more than two 
degrees of freedom are involved. The Sylvester method 
of obtaining the condition that two simultaneous equa- 
tions have a common root completely eliminates the 
unknown quantity. It is feasible, however, to termi- 
nate the process of elimination with two equations of 
the first or second degree. The choice made in the 
following sections is the use of two equations of the 
first degree. 

The equations arising in the calculations in the case 
of three degrees of freedom are of the form: 



A 3 X 3 +A 2 X 2 +A 1 X+A Q =0\ 
B z X*+B 2 X 2 +B l X+ B 0 =0\ 



(1) 



where in special cases the degrees of the equations [(3,3) 
in equation (1)] may be (3,2), (2,2), (2,1), or (1,1). The 
quantity X is an unknown frequency parameter, and 
the coefficients A and B are functions of a large number 
of parameters: structural parameters a, J, c, x ai x^, r a 2 , 
rf* *> 9 a i 9& and 9»t frequency parameters ft*, Q«, Q fi ; 
and the reduced frequency Ifk. For a particular air- 
craft structure represented by given parameters there 
corresponds a flutter velocity and a frequency deter- 
mined from X and l/k. Expressions for the quantities 
A and B are listed in references 1 and 2. In the follow- 
ing discussion it is assumed that these quantities are 
available. 
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The common solution of equations (1) can be ob- 
tained from the common solution of 



brX+b { 



Oo=0\ 



(2) 



where a u do, b u and b 0 are functions, listed later, of the 
^4's and B's in equations (1). Now, from equations (2) 
it is evident that the common solution exists if and only 
if 

X x ——a^/a l is also equal to X 2 = z —b Q /b l 

Then, if all the parameters but one are kept constant, 
for instance l/k, and X x and X 2 are plotted against 1/k, 
the intersection (or intersections) determines the com- 
mon root (or roots) X and the value (or values) of 1 fk 
for which this common solution occurs, and A" and l/k 
together determine the flutter solution for the particu- 
lar structure. 

Another possibility, namely, keeping l/k fixed and 
plotting X against one of the structural or frequency 
parameters, will yield as a flutter solution the necessary 
structural parameter. Many variations are possible. 

The Sylvester resultant of equations (2) is the deter- 



minant 



b\ b 0 



and its vanishing is the condition for 



the existence of a common root. If this quantity is 
plotted against l/k as the abscissa, for instance, the 
intersection with the 1 \k axis gives the required value 
of l/k. The first-mentioned method involving two 
parameters is preferable, however, because the two 
curves are simpler and yield both X and l/k simul- 
taneously. 

There remains, then, only the task of listing the 
expressions for a u a Qy b u and b 0 . It is convenient to 
list these expressions separately for the cases in which 
the degree of the equation is (2, 2), (3, 2), and (3, 3). 

In order to obtain do, a u b Q , and 6 X for the case of two 
quadratics, multiply the first of equations (1) by B 2 
and the second equation by A 2 and subtract; and 
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similarly multiply the first of equations (1) by B : X+Bx 
and the second by A s X+Ax and subtract. Then 



a 0 = 



A 
Bo 

Ax 
Bx 

A 0 
B 0 

A 0 
Bo 



A 2 
B 2 

A, 
B 3 

A, 
Bx 

A, I 
B 2 1 



= a 0 



Similarly, for one cubic and one quadratic (3, 2) : 

a 0 = 



0,1- 



B 0 



60= 

l 

In the case of two cubics 



a 0 = 



AoBi-AiBt -BoA 3 
B 0 B 2 

B t B 2 

A Q B l -A l B 0 AoBi-AtB* 



Bi 



ai= 



60= 



bi = a 0 



A 0 


A, 


\A 0 


A 3 


Bo 


Bx 


i 5 0 


B 3 


A 0 


A, 


! A 2 


A 3 


B 0 


B 3 


i B 2 


B 3 


A 0 


A 2 




Ao 


A 3 


Bo 


B 2 




Bo 


B 3 


A x 


A 3 




A 2 


A 3 


Bx 


B 3 




B 2 


B 3 


A 0 


A t 




Ao 


A 2 


Bo 


Bx 




Bo 


B 2 


Ao 


A 3 




A x 


A 3 


Bo 


B 3 




B x 


B 3 



In the use of this method it is sometimes found that 
the common intersection is not obtained with precision 
without the use of many values of \jk. It may then 
appear to be more convenient to employ a different 
form. Thus, in the case of two cubics, there are three 
possible forms for a 0} a u b Q , and b x and a second form is 



a Q ~ 



b 0 = 



A 0 


A 3 


B 0 


B 3 


Ao 


A 2 


Bo 


B t 




Ax 




B l 


Ao 


A 3 


Bo 


B 3 


Ao 


A 3 


Bo 


B 3 


Ao 


A t 


Bo 


B t 



A 2 A$ 

B 2 B$ 

Ai A z 

Bi B z 



+ 



A t A 2 
Bi B 2 



A 2 


A 3 


B t 


B 3 


Ax 


A } 


Bx 


B 3 



A x A 3 
Bx B 3 



I Ao 


A 3 




Ax 


A 2 






+ 






B 0 


B 3 




Bx 


B 2 



The method is not limited to the original form of the 
equations. Assume, for example, that both X and l/k 
are preassigned and that it is required to know the 
values of two parameters, say u x and t^, which have as 
the flutter solution the preassigned values of X and 
1/k. The original equations can be considered as 
equations in u x and u 2 whose common solution is 
determined by 

aitti+a 0 =0 
Mi +60=0 

where a u a 0f b u and b 0 are known (calculable) functions 
of all the other parameters. If the two roots are 
plotted against W2, the intersections (if any) will give 
the required values of i/i and v*. 



APPENDIX C 



EXAMPLE ON FLUTTER OF BIPLANE 



Experiments on the vibration frequencies showed the 
following results (values given in cycles per min): 

1. Antisymmetrical torsion of wing-cellule 

system 1300 

2. Symmetrical bending of wing-cellule sys- 

tem 800 

3. Symmetrical torsion of wing-cellule system- 1300 

4. Local wing bending: 

a. Lower wing, with node at or near inter- 

plane strut 1300 

b. Upper wing, with node at or near inter- 

plane strut 1100 

5. Aileron against controls 1100 

6. Local torsion in aileron 1S00 

7. Local torsion in flap 1100 

8. Engine rocking 830 

There are two possible types of ternary flutter: 

a. Symmetrical torsion-symmetrical bending- 
symmetrical aileron motion. The frequencies are 1300, 
800, and 1100, respectively. 

b. Antisymmetrical torsion-antisymmetrical bending- 
antisymmetrical aileron motion. The frequencies are 
1300, 0, and 0, respectively. 

The other parameters were used as follows: 
a—— 0.2 (elastic axis at 40-percent chord); x a ~0.2 
(center of gravity at 50-percent chord; the actual center 
of gravity was near 48-percent chord); r a 2 =l; *=0.2 
(this value of the wing-density parameter corresponds 
not to sea level but to an altitude of approximately 
10,000 ft); 2^=4 feet 9 inches (reference chord). 

With the use of these parameters, there is obtained 
for the torsion-bending (case 1) flutter-speed coefficient 
v/bo) a from figure 1 a value of 1.26. The reference 
velocity bo> a is equal to 221 miles per hour. Thus the 
flutter speed V F is equal to 278 miles per hour. Because 
the observed flutter speed on this biplane was lower 
than this value, (about 200 mph), the aileron was evi- 
dently involved. The parameters relating to the aileron 
were assumed to be as follows: 

Location of the center of gravity, xp 0. 002 

Radius of gyration, rf 0. 002 

Chord location, c _.. 0.6 

The aileron was considered a full-span aileron. This 
assumption is fairly reasonable because the lower wing 
flap was almost identical with the aileron. These 
values were used in the results shown in figure 1, which 
16 



was based on the biplane. The ratio =0.833 
gives, for the assumed unbalance r$ — 0.002, a value of 
the flutter coefficient Vjbw a of 0.68 or a speed of 151 
miles per hour. 

For the antisymmetrical case, if a full-span aileron 
and zero damping are conservatively considered, there 
is obtained from figure 37 the value r/6w a =0.41. A 
value of the internal damping g a of 0.01, however, in- 
creases the flutter coefficient, to 1.18, which is equal 
to 261 miles per hour (true speed). Notice that this 
value is calculated without the benefit of a fractional 
aileron. If there is used in the symmetrical case a 
small value of the internal damping g a of 0.01, it is 
seen from figure 2 (b) that there is only a slight favor- 
able effect from mass-balancing: The flutter coeffi- 
cient v/beca is equal to 1.10 for x 3 ~ 0.002 and increases 
to 1.16 for x $ = — 0.002. With the use of r/iu« = l.l, 
there is obtained a flutter speed of 243 miles per hour 
(true speed). From later experiments it has been 
found that the value g a = 0.0l is evidenth" a safe value 
to use in such calculations. It is thus noted that the 
flutter speed, because of this effect, approaches the 
torsion-bending value. It is further observed that with 
this amount or a larger amount of damping the mass- 
balancing of the aileron becomes fairly ineffective. 

Since the calculation for the symmetrical case based 
on 0 a = O.Ol gives values of the flutter velocity in the 
order of 240 miles per hour, true speed (corresponding 
to an indicated speed of approximately 206 mph), it 
is probable that this case describes the observed flutter, 
which was known to be symmetrical. 

This biplane was aerodynamically cleaner than many 
of the earlier types and it is possible that the absence 
of numerous interplane wires and struts contributed to 
a lowering of the torsional damping effect to such an 
extent that flutter was invited. No doubt, many of 
the older types of biplane were safe from flutter because 
of their large structural damping. 
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TABLE I 
[tJ, I; k, 0.2; c, 0.6; a, -0.2; jr^, 0.2; r i, 0.002; £, 1.0) 
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TABLE III 
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.607 
.607 
.607 

.007 

.607 



0.2 



0.5 
.8 
.9 

1.0 

Variable 
.8 

1.0 
.8 

1.0 
.8 

1.0 

1.0 



.01 
.01 

0 
0 

0 

.01 



TABLE IV 

[rj, 0.5; «, 0.2; c, 0.6; a, -0.2; x m , 0.2; rj, 0.002; {, 1.0] 



r igurc 







J > 












0. 60 1 


— 0. 005 


0 


0 


0 


12 


do 


. 607 


— . 002 


0 


0 


0 






. 607 


0 


0 


0 


0 




'""do""!"."'!" 


. 607 


.002 


0 


0 


0 


13(a)... 


/....do 


.607 


.002 


0 


0 


0 


13(b) .... 


1 - do 


.607 


.002 


.01 


0 


0 


/--■do 


.607 


0 


0 


0 


0 




t ill*. . . 


. 607 


0 


.01 


0 


0 


13(c)--. 


/....do 


.007 


-.002 


0 


0 


0 




1 ...do 


.607 


-.002 


.01 


0 


0 


13(d).... 


/....do 


.607 


-.005 


0 


0 


0 


14 (a)... 


\ do- 


. 607 


-. 005 


.01 


0 


0 


0.607 


.607 


. 002 


Variable 


Variable 


Variable 


14(b)... 


.833 _ 


.607 


. 002 


Variable 


Variable 


Variable 


15(a).... 


/Variable 


.316 


.002 


0 


0 


0 




\....do 


.316 


.002 


.03 


0 


0 


15(b)....| 


/----do 


1.0 


.002 


0 


0 


0 


15 (c).... 


1. ...do. 


1.0 


.002 


.03 


0 


0 


/....do. 


. 316 


-.002 


0 


0 


0 




1 do 


.316 


-.002 


.03 


0 


0 


15(d)....i 


(....do 


1.0 


-. 002 


0 


0 


0 


i 


U--do..__ 


1.0 


-. 002 


.03 


0 


0 



TABLE V 
[rj, 0.5; c, 0.6; a, -0.2; rf, 0.002; |, 1.0J 
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TABLE VI 
IrJ, 0.25; c, 0.6; a, -0.2; r m , 0.2; r^, 0.002] 



Figure 


P ■ 






ff 


( 




- 

29 


/Variable 


0. 607 


0.2 


0.002 




1.0 


— — 

0 


»....do . 


.607 


. 2 


-.002 


1.0 


o 


30(a) 


/....do 


.607 


. 2 


.002 


I. 0 


o 


{....do 


.607 


. 2 


.002 


1.0 


.01 


30 (b) 


/... do 


.607 


. 2 


-.002 


1.0 


o 


\....do 


.607 


. 2 


-.002 


1.0 


.01 




/....do 


.607 


.2 


.002 


1.0 


() 


t-.-.do . . . 


. 607 


. 2 


.002 


.8 


o 


32 (8) 


f... do . 

{....do.;.::... :: 


. 316 
. 316 


. 2 
. 2 


.002 
.002 


1.0 
1.0 


o 

. 03 


32 (b) 


j do 
1 <lo 


I. 0 
1. 0 


m 2 
m 2 


. 002 
. 002 


1.0 
1.0 


o 

03 


32(c) 


)'..Ao..'.'. " 


!316 


*2 


-!o02 


1*0 


0* 


\....do 


.316 


.2 


-.002 


1.0 


.03 


32 (d) 


/....do 


1.0 


.2 


—.0(12 


1.0 


0 


\....do 


1.0 


.2 


-.002 


1.0 


.03 


33 (a) 


/....do 


.316 


. 125 


.002 


1.0 


0 


\ ...do 


. 316 


. 125 


.002 


1.0 


.03 1 


33(b) 


/....do 


.607 


. 125 


.002 


1.0 


0 


\---do 


.607 


. 125 


.002 


1.0 


.02 I 




(.-.do 


1.0 


. 125 


.002 


1.0 


o ! 


33(c) 


|....do 


1.0 


.125 


.002 


1.0 


.05 | 




I... .do. 


1.0 


. 125 


.002 


1.0 


.03 


34(a) 


/....do... 


.316 


.125 


-.002 


1.0 


0 


\....do 


.318 


.125 


-.002 


1.0 


.03 


34 (b) 


/....do 


.607 


. 125 


-.002 


1.0 


0 


\....do 


.607 


. 125 


-.002 


1.0 


.01 


34 (c) 


/....do 


1.0 


. 125 


-.002 


1.0 


0 


1 — do... 


1.0 


.125 


-.002 


1.0 


.03 



TABLE VII 

{rj, 0.25; *. 0.0752; c, 0.5; a, - 0.4; r., 0.2; x fi , 0; r 0 \ 0.002] 



Figure 



35(a). 
35(b). 
35(c). 
35(d). 



Variable . 

do... 



/.-..do.. 
\ ...do., 
/...-do.. 
l....do . 

;/__. do 
:\...do.. 



0. 316 
. 316 
. 500 
.500 
. 707 
.707 

1.0 

1.0 



TABLE VIII 

\rj. 0.25; c, O.G; o,- 0.4; x„. 0.2; x fl , 0; r fi t, 0.0012] 



36 (a)... 
36(b)... 



/Variable 

\.--.do 

|-.--do 

....do 

l.-.-do 



0.25 
.25 
.25 
.25 
.25 



0.25 
.25 
. 1 
. 1 
.1 



P. 


n 


<\ 


0 


0 


0 


. 30 


. 10 


. 10 


0 


0 


0 


.02 


.02 


.02 


. 10 


. 10 


. 10 



TABLE IX 

fws/«„ 0; « k '« al 0; «, 0.2; c, 0.6; a, - 0.2; r m , 0.2; iy, 0.002] 



j Figure 






* 


fa 






0.004 

.002 


1.0 

1.0- 


Variable. 
Do. 






-.002 


1.0 


Do. 
Do. 

. 0 






-.006 

.004 


1.0 _. 

Variable. 


38 




.002 


...do 


0 






-.002 

-.006 


....do. 

....do 


0 

f 0 






Variablo 


.5 


1 0 






...do 


.6 


0 






...do 


.7 




39 




....do 


.8 


1 0 

! 0 






.-..do.. 


.85 






...do 

....do 


.9 

10 


0 

0 




' .5 


.002 


1.0... 


Variable 
Do. 
Do. 




.5 


-0.002 


1.0 


40 




.002.. 

-0.002 

.002 

-.002 


1.0 


I 


.25 


1.0 

1.0... 

1.0 


Do. 
Do. 
Do. 
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